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We compute the drag force exerted on a quark and a di-quark systems in a background 
dual to large-N QCD at finite temperature. We find that appears a drag force in the former 
setup with flow of energy proportional to the mass of the quark while in the latter there 
is no dragging as in other studies. We also review the screening length. 
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1. Motivation and results 



QCD at T 7^ plays an important role in understanding two related areas in physics: 
the physics of the early universe and the physics of heavy ions collisions. Due to the 
increasing experimental evidence for the existence of the quark-gluon plasma (QGP) state 
at RHIC and the existing planned experiments as LHC and FAIR we shall be 

concerned with the latter. The experimental results are against the old prejudices that 
assumed that at sufficient high temperatures the interactions between quarks and gluons 
would lie in the perturbative regime making of the QGP an almost ideal gas. The data 
suggest that in a QGP matter is almost, but totally, in the deconfined phase but with the 't 
Hooft coupling constant in the ballpark of a decade invalidating any perturbative approach, 
Qym^ ~ 10. This makes of QGP a good testing ground of the AdS/CFT duality that has 
as a prerequisite the condition Qym-^ 1- particular, to model this physics, we shall 
use a setup of D4 branes at finite temperature [|5| but with a zero chemical potential. 
This, in principle, does not present an obstacle to attain realistic results inasmuch we are 
just interested in the fluid properties produced in the central rapidity region 0. 

We shall be concerned mainly with one measured quantity, the energy lost of a high- 
energy parton in the QGP 0. We look at two different systems: a single heavy-quark and a 
di-quark passing through the QGP. To make connection with string theory one models the 
single external quark as a string dangling from the boundary and the di-quark system as a 
string hanging between the two quarks. Both setups have been extensively studied recently 



in the same context but with different backgrounds [l8|,|9|,p!0i 
These strings extend on the radial transverse coordinate and can be thought as a flux tube 
spreading out of the 3-1-1 boundary theory. This flux tube will travel along the bulk that 
is constituted by a thermal bath of gluons. The interaction with this gas will exert a force 
on the string that is propagated along it to the external quarks at the boundary. This is 
the effect we are precisely interested in. 

As we shall explicitly see there are no qualitative changes in our results to what is 
already known in the literature: a single quark experiences a drag force while in the di- 
quark system there is no dragging in the direction of the movement. These seem to be the 
general pattern to all the models studied so far. 

The paper is organized as follows: sec. 2 contains the two supergravity backgrounds 
dual to large-N QCD we are considering, one at zero temperature and the equivalent at 
finite temperature. In sec. 3 we study, in both phases, the force experienced by a heavy- 
quark passing through a plasma and discuss the thermal mass. Sec. 4 contains the study, 
in the decofined regime, of a pair of quarks. We focus the attention in the energy-distance 
function. Because is not clear so far whether from this or similar studies can be obtained 
the jet-quenching parameter itself, q, we retreat to go further and obtain it from our 
results. Furthermore an analysis on the stability |j21[] or favored energetic configuration 
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2^ is lacking in our study. 



2. Zero vs. finite temperature QCD at large-N 

We briefly review the two backgrounds we used. In doing so we must bear in mind 
that the field theory we expect to describe is at the boundary of both backgrounds and 
therefore they must asymptote one to each other in this limit. The reason to use two 
different setups is related to the confined- deconfined phase transition. At zero chemical 
potential, we expect a cross-over region between these two phases in which hadrons behave 
so differently. Notice, that we do not use probe-branes in our analysis and hence, we are 
mainly testing properties of the thermal bath while matter is treated as a external source. 

At zero temperature we use a stack of D4-branes (both models are in Lorentian sig- 
nature) were the initial time coordinate has been compactified on a circle of radii u^- This 
breaks supersymmetry and give masses to the fermions at one loop. After this we Wick 
rotate one of the remaining space coordinates that will play the role of time. The final 



form of the various fields: metric, the RR four- form and the dilaton, are \22> 
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Notice that as the temporal coordinate can not develop a singularity after compactification 
the horizon coordinate is not related with the temperature of the field theory, and we are 
free to choose T = 0. The peculiarities of (|2.1| ) are already well studied in the literature: 
i) It confines, ii) It has chiral symmetry breaking when probe-branes are added, iii) There 
are indications for the decoupling of the Kaluza-Klein modes from the spectrum. 

The second of our backgrounds is an euclidean continuation of ( |2.1| ) and is assumed 
to describe physics above the critical region, T > Tf.: time is compactified over a circle 
of period (3. Anti-periodic boundary conditions for the fermionic content of the theory is 
imposed along t and X4-directions. The latter with radii 2t{R. 

The field content is 
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where to avoid time direction singularities one must choose 
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Because there is no smooth interpolation between the two backgrounds, ( |2.1D -(PT^), 
the system is believed to describe a first order transition. By looking at the free energy 
one sees that (|2.1| ) is the relevant background at zero temperature while (|2.2| ) is at high 
ones P,^. The relevant salient point of ( |2.2| ) is that while it supports a discrete hadron 
spectrum, chiral symmetry breaking is still at work, it has occurred deconfinement, the 
gluonic sector has lost its symmetry. Thus all in all, ( p.2| ) can contain quarks, gluons and 
hadrons simultaneously as degrees of freedom. This picture, although not very intuitive, 
is in agreement with the QGP data where at T = not all the hadronic dof are lost [H . 



3. Drag force for a single quark 

We focus first in the study of the force exerted by an external agent on a quark 
1^,^. For that purpose we fixed all but the five uncompact dimensions, from which all 
transverse coordinates vanish except one x. Furthermore we work on the static gauge 
t = T,u = a, x(r, a). The action of a fundamental string in the string frame is given by 
the Nambu-Goto action 

If 

S = —- — / dadr^/^ . (3.1) 



2TTa _ 

3.1. Deconfined phase 

With the previous, static gauge, ansatz the lagrangian density for ( p^) reads 



where dots (primes) stand for derivatives w.r.t. r (a) . We have defined h = (^) and 
dropped all the arguments in the functions. The e.o.m. for x is 

a(^)-fajn;^)=o. (3.3) 



-9J \V-9. 

Notice that the solution of a static string, x = ct, stretching from the horizon up to the 
boundary exists. By choosing x(r, a) = vt + ^(cr) in ( |3.3| ) one obtains in terms of the 
worldsheet conserved quantity, tt^, 

(3.4) 
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This expression is cross-checked with the generic, diagonal metric, solution of The 



expression (|3.4| ) is much like the equivalent one in A/" = 4 SYM . Demanding that any 



piece of the string has a subluminal velocity amounts to fix the sign of the squared root in 
( p.4|) at certain critical point u^^ = Uj-/{1 — v'^) . This leads to the fixation of the constant 
of motion 

.,^r-L-f''^. (3,5) 



V 



The main feature of the latter expression w.r.t. the A/" = 4 case is the equal behaviour in 
terms of the velocity. This probably signals an universal type of behaviour, independently 
whether the model is derived from a supersymmetric one or not. Inserting ( pi .51 ) in ( |3.4D and 
using ( |2.2| ) leads to a remarkably simple relation 

The shape of ( |3.6| ) in the x — u plane at a given time is a monotonic bended curve. One 
could ask whether this bending signals the interaction of the quark with the reservoir of 
gluons. To answer this we look at the density of energy and the a;-component momentum, 
( p.2|) does not depend explicitly neither on t nor in x, obtained from the conserved currents. 




l_ f l + hf{x'Y \ nl\ _J_ vfhx' 



(3.7) 



Inserting the configuration ( p.6|) in (^]3) we find that, as in the A/" = 4 case, the two 
currents tt^, tx] are constant along the string 
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The features of ( |3.8| ) are: i) it vanishes for f ^ indicating that no force is needed if 
the quark is kept at rest w.r.t. the fiow. ii) If we pull the quark with constant velocity, 
transmitting to it energy at the ratio 7r|, this energy goes along the string were it is 
dissipated into the media. The fraction of energy supplied to the surrounded thermal bath 
at a given point along the string is constant. If one calculates the total of energy supply 
to the system, danj:, turns to be infinite. This is clearly understood if we bear in 
mind that we deal with infinitely massive quarks, and hence an infinite amount of energy 
is needed in order to move it. As a matter of fact, without the addition of probe-branes 
|25| j to hold the quark at finite distance from the horizon, the only stable solution is to 
set the quark at infinity. If we were dealing with the latter set up, finite quark masses, the 
lost of energy would be finite at any point, and the total energy needed to pull the quark 
would be proportional to its mass (location along the radial coordinate). 
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3.2. Confined phase 

We comment briefly on the setup at zero temperature. The expectations are the 
foUowing: at zero temperature and a strong coupling constant we know that quarks are 
conflned, they are presented in singlets of color. It is mandatory then, that an embedding 
as the previous one, a single quark with velocity v, fails to be stable. The lagrangian 
density in the static gauge reads 



1 liu'y 



27ra V / 
that leads to the following e.o.m. 
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As in the parallel previous case, (|3.9|) does not depend explicitly on x and there is a 
conserved worldsheet current, tt^. Solving ( p.lOj) in terms of these we have 
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(3.11) 



Notice that in the limit / ^ 1 (|3.9| )- (|3.11J ) reduce to the equivalent expressions of the 
previous section as it should be. 

Let's inspect ( p. 11 ). The functions / and h are strictly positive in their domains. 
Hence there is no way of flxing the conserved quantity tc^ as was done in ( p.4|) : the critical 

2/3 

value for the radial coordinate is Ucr = Rtc^ , thus as we increase the momenta the 
penetration of the string in bulk also increases, but at finite Ucr the velocity of the string 
turns to be supraluminical. This picture is consistent signaling an unstable configuration. 

Just for ending this subsection, and it would be proven useful below, we quote the 
conserved currents 
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that also fulfills the relation tt/ = —vivl 



vhx' 
—hx' 



(3.12) 



4. Thermal mass 

We compute the energy (mass) for a static configuration in ( p.2|) and check whether 
it satisfies the field theory expectations as a function of the temperature. 
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Using (|2.1|) and the embedding t = t,u = a,x{a,T) = xq the mass at zero temperature 
is given by 

m{T = 0)=E = - ^dan^ = ~Rl^+Urr.2F^(-ill-^^) , (4.1) 



9 r(i) ' "'-^ ^ V 3'2'3'729z.3^/ 

with Um a hard cut-off that regularizes the expression. 
For (|2.2| ) one has a very simple expression 

m(T) =E = - dan^ = ~ -^-f- . (4.2) 

Both expressions, ( [1.1|) and ( [4.2| ), differ only in a constant in the asymptotic regime. 
This is due to the difference in integration limits. As mentioned above the only stable 
solution without introducing probe branes is to locate the quark at the boundary, Um oo. 
This makes an infinitely heavy-quark. A possible way to renormalize (|4.2| ) is to substract 
( [4.1|) . This will lead to the value of the thermal mass at the boundary. The general pattern 
of the mass as a function of the temperature is as follows: For low-temperature ( [4.1| ) is 
bigger than (|3). At /3 = 2TT^/^R^yT (1/6) /r(2/3) both become equal and for higher 
temperatures they exchange their roles. Notice that in increasing the temperature, up to 
its maximum value (|2.2D the mass, m(T) — m(T = 0), decreases in accordance with the 



FT expectations [p6| , |27[| . This corroborates in an independent way some of the findings in 

m. 



5. Drag force for Quark- Antiquark system 



In this section we shall only consider a di-quark system at high temperatures. The 
configuration we shall focus on is based in the studies of • The parallel of for 

the background (|2.2|) has been studied thoroughly in . In the static gauge we shall use 
the ansatz t = r , x{a, r) = vr + ^(cr) , y{a) , z{a, r) = a , i.e. the velocity is perpendicular 
to the dipole moment. More general ansatz can be studied ||2^,0 but in view of our results 
we do not expect qualitative changes w.r.t. them. 

First we consider some general features of ( p.2|) w.r.t. this embedding ||3^. The 
lagrangian density reads 



-7 



l + fh{y'f +fh{a 
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(5.1) 



From the e.o.m. we obtain the trajectory in terms of worldsheet associated currents 
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f^hif-v^){fh-nl)-fhnl 



1 



hif-v^){fh 



(5.2) 



7 



As in (|3.4D , due to the signature of the factor inside the squared root, exists a condition on 
the minimal value of the radial coordinate. In the case of a string hanging from a quark 
was u > UcT- Using (|5.2|) we find that Wmin > signaling that the string stretching 



between quarks can not penetrate up to the horizon. The equality, Wmin = Ucr , only 
holds if TCy = 0. As a last remark we noticed that the only possible string configuration 
has the string parallel to the ^-axis: the condition for a derivable curve at rmin imposes 
dx / dy\r^.^ — > oo that only is posible if tTj; — > 0. As a consequence, x' = and the string 
moves with constant velocity and with a shape parallel to the 2;-axis. 

After these remarks we restrict the above ansatz to the case ^(a) = 0. In order to 
obtain the effect of a possible drag force one could boost the system from the plasma 
reference frame to the di-quark system one: x ^i-l^ ' ^ ~^ Vi-v'^ ' -^^^^ 
customary to do so one should bear in mind that the experiment is in the plasma frame, 
and we stick to it in the remainder. Whether we perform the boost or not the profile of a 
string with velocity v along the a;-axis and extended in the radial coordinate z is unaltered 
because the boost is perpendicular to its shape 

y' = ±R\^ ul.^ ^ T y . (5.3) 



The end points separation of the strings leads at a distance L one from each other 



L{ny,v)=2 duy'. (5.4) 



To perform the integrations we follow the prescription in |]3^ and change variables to 
h = 1 — (ut/u)^ . In addition we also make use of the rescaled momenta fy = iXy [R/ut)^ . 
With these ingredients the integration domain becomes [{fy + f ^) /(I + ), 1]. After some 
algebra, and using (|2.2| ), (|5.4|) reads 



L{fy,v) = ^^ dh . (5.5) 

27rT (1 _ hy/^^h {f^{h -l) + h-v^) 

This integral can be casted in terms of hypergeometric function 

r(f) ^ A 1 7 i + f^ 



Despite its appearance its domain remains in the reals. 



P2 
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Fig. 1: (r.h.s) Diquark system separation in unit of l/(27rr) . Prom top to bottom 
(blue) curves displays the values v = 0,0.45,0.7 and 0.99. The (green) thicker 
vertical curve is the fit to the function I/max(/y, ■«). (Ihs) depicts, with the data in 
boxes, the maximum distance as a function of the velocity. The green curve is a 
quadratic fit. The triangles are the screening length. 



In the l.h.s of fig. 1 we have plotted the quantity / = 2nTL as a function of the 
external force fy. The results are much like the same as in the supersymmetric case 

and also is in accordance with the static one the separation between quarks 

is not a monotonic function of the external force. As a consequence for a given velocity 
V there is a maximum separation, Lm_a.x{fy,v). This function can be written analytically 
in terms of hypergeometric functions, but its actual expression is not very illuminating. 
Instead of presenting it we have added to the figure a few representative points, and a cubic 
polynomial fit to them. In the r.h.s of fig. 1 we have evaluated the maximum separation 
between quarks as a function of the velocity. Points belong to the same set as in the r.h.s 
figure and the green curve is a fit. The best x-squared is obtained with the function 
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For V > 0.98 there is a pronounced drop in the curve and it vanishes for v = 1. 

We turn now to the evaluation of the string energy in the plasma rest frame. Using 
the density lagragean (^TT]) reduces to 



(1 



Un 
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(5.7) 



and the string action reads 

S {TXy.V) 



1 



na'T 



r/2 



r/2 



dtdr C 



(5.8) 



The above expression, ( ^.8|) , diverges and needs to be regularized. This is more easily seen 
by looking the behaviour of (^.7|) at the boundary: Vl — f^, which after the worldsheet 



9 



integration per boundary unit time of ( ^.Sp becomes linear divergent. One possibility to 
cure this is to subtract the action of a pair of quarks stretching from the boundary up to 
the horizon 



1 



S{'Ky,v) = — — / dt\\ drC- drC] . (5.9) 



T/2 \J Umin ''Ut 



This quantity is related with the expectation value of the Wilson loop operator. Even if it 
is an interesting quantity at zero temperature, in pure gauge theory at finite temperature 
is the Polyakov-loop the order parameter which signals whether the system is confined or 
not For this reason instead of evaluating ( pT9| ) we compute the energy stored in the 
di- quark system. 

Using (|5.8| ) one can already compute its associated conserved canonical momentum 
densities 

As one can see the component tt^ vanishes signaling that, contrary to the single quark 
configuration, the di-quark experiences no force by the presence of the gluon bath in the 
direction on its movement, in that sense one can recover the experimental data that signals 
the dynamics of quarks inside an almost perfect fiuid. Notice eventhought that exists a 
force in the y-direction.This behaviour is reminiscence to the image of the Cooper pairs in 
super-conductivity: one electron can feel a drift while a pair does not. 

As already happened with ( |5.8| ) the energy of the system will need to be renormalized. 
We shall chose the standard method subtracting the bare mass of a pair of quarks, but 
other scheme, as for instance subtracting the corresponding energy obtained from ( p.l|) , as 
was done in sec. 3, will do the same job as far as it has the same asymptotic behaviour. 
The final result reads 

(l-/l)4/3//2(/, ^2 Jo (1-/1)4/3^/1^ 

(5.11) 

In the asymptotic regime this expression vanishes. More over due to an infrared singularity 
( ^.11|) is not uniquely determined ||3^ . 



We have plotted in the l.h.s. of fig. 2 a normalized energy, e = 3E/{ut), vs. the 
separation of a quark pair at different velocities, that turns to be a bivaluated function. 
Bound states are allowed to exist only in the region were e < 0. If e > the energy 
of a pair of free quarks is less than the energy of the bound state and chiral symmetry 
restoration has occurred. The separation between quarks where this phenomenon happens 
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Fig. 2: (r.h.s.) Energy as a function of the quark separation for different velocities. 
From top to bottom v = 0.8, 0.7, 0.4 and 0. (l.h.s.) The maximum energy ( |5.11 ) as 
a function of the velocity. Data points are take from the r.h.s. plot. The curve is 
a quadratic fit to the point with v < 0.6. 



is defined as the screening length (L^) as increasing the distance between quarks 

gluons start to screen the quark interaction until it becomes negligible and deconfinement 
has occurred. Notice that if the system has some velocity v the screening length diminishes 
very smoothly. In the static case one finds ~ 1.34/ (27rT), slightly lower that the J\f — 4 
case. The relevant feature of fig. 1 r.h.s. is that the screening length is almost a constant 
factor as a function of the velocity for f < 0.9 in agreement with the conclusions in [^^| 
where it was spelled out that the main velocity dependence comes from the boost action. 
In this region never reaches the maximum possible separation between quarks, i.e. 
matter deconfinement occurs at shorter distances than naively expected, a factor ~ 3/2 
lower than the expected. From v > 0.9 on, v 1, the screening length, L^, is given by 
the maximum separation between quarks and goes to zero. 

Inspecting more carefully fig. 2 l.h.s. the salient points are: i) For each separation 
between quarks we find two possible energetic states: a lower energy position, that is the 
stable one, and an upper that indicates some metastable state. For v >> and for high 
energies the curve is almost fiat with a very large energetic gap between the two possible 
states given a di-quark separation. This indicates that the metastable states are very 
unstable. For all the curves, the rhs part of the plot, very short separations, can only 
be reached for — > cxd. ii) For a fixed energy, there are also two possible positions, the 
one with the minimum separation between quarks correspond to the set of metastable 
states. Although not very intuitive quarks do not lie as near as they can from each 
other without decaying, in) As we increase the velocity the curve moves to the right and 
upwards. One should wonder about the different shape and velocity behaviour for the cases 
of medium-high velocities in fig. 2 w.r.t. the equivalent ones in [^. We have checked that 
this difference comes from the boost, if we go to the di-quark rest frame the qualitative 
behaviour between the two models agree. 

In the r.h.s of fig. 2 we have plotted the maximum energy of the di-quark system at 
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a given velocity. Up to v ~ 0.9 this maximum energy is reach at the maximum separation 
point between quarks and is an increasing function. At this precise point, v ~ 0.9, the 
maximum of energy is not longer found at this distance and we find a decreasing function 
that diverges at v — > 1. 

6. Summary 

We have explored the drag force for a single and a di-quark system in a dual model 
to large-N QCD. Particularly relevant is that the gross features we find in this model 
are in qualitative agreement with those derived in others models obtained from A/" = 
4. This signals that the main physical aspects are not sensible to the initial model we 
depart at the very beginning, before breaking supersymmetry incorporating temperature, 
is super symmetric or not. 

Our results for the single quark are parallel to previous studies. The behaviour of 
the thermal mass as a function of the temperature is in agreement with the field theory 
expectations. 

In addition the di-quark system shares most of the peculiarities of other systems 
already studied in the literature. The main relevant difi'erence coming from the almost 
constancy, up to velocities v < 0.9, of the screening length as a function of the velocity 
and its relative smallness in from of the maximum allowed distance between quarks. 
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